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Abstract 

One of the most remarkable recent developments in the study of ultracold 
Bose gases is the observation of a reversible transition from a Bose Einstein 
condensate to a state composed of localized atoms as the strength of a periodic, 
optical trapping potential is varied. In [I] a model of this phenomenon has been 
analyzed rigorously. The gas is a hard core lattice gas and the optical lattice 
is modeled by a periodic potential of strength A. For small A and temperature 
Bose-Einstein condensation (BEC) is proved to occur, while at large A BEC 
disappears, even in the ground state, which is a Mott-insulator state with a 
characteristic gap. The inter-particle interaction is essential for this effect. This 
contribution gives a pedagogical survey of these results. 
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1 Introduction 

One of the most remarkable recent developments in the study of ultracold Bose gases 
is the observation of a reversible transition from a Bose-Einstein condensate to a state 
composed of localized atoms as the strength of a periodic, optical trapping poten- 
tial is varied [21 Oj- This is an example of a quantum phase transition j3] where 
quantum fluctuations and correlations rather than energy-entropy competition is the 
driving force and its theoretical understanding is quite challenging. The model usu- 
ally considered for describing this phenomenon is the Bose-Hubbard model and the 
transition is interpreted as a transition between a superfluid and a Mott insulator 
that was studied in [H] with an application to He 4 in porous media in mind. The 
possibility of applying this scheme to gases of alkali atoms in optical traps was first 
realized in j^]. The article |7j reviews these developments and many recent papers, 
e.g., jHISMmilSII^dllSlini are devoted to this topic. These papers contain 
also further references to earlier work along these lines. 

The investigations of the phase transition in the Bose-Hubbard model are mostly 
based on variational or numerical methods and the signal of the phase transition is 
usually taken to be that an ansatz with a sharp particle number at each lattice site 
leads to a lower energy than a delocalized Bogoliubov state. On the other hand, there 
exists no rigorous proof, so far, that the true ground state of the model has off-diagonal 
long range order at one end of the parameter regime that disappears at the other end. 
In this contribution, which is based on the paper pQ, we study a slightly different 
model where just this phenomenon can be rigorously proved and which, at the same 
time, captures the salient features of the experimental situation. 

Physically, we are dealing with a trapped Bose gas with short range interaction. 
The model we discuss, however, is not a continuum model but rather a lattice gas, 
i.e., the particles are confined to move on a d- dimensional, hypercubic lattice and 
the kinetic energy is given by the discrete Laplacian. Moreover, when discusssing 
BEC, it is convenient not to fix the particle number but to work in a grand-canonical 
ensemble. The chemical potential is fixed in such a way that the average particle 
number equals half the number of lattice sites, i.e., we consider half filling. (This 
restriction is dictated by our method of proof.) The optical lattice is modeled by a 
periodic, one-body potential. In experiments the gas is enclosed in an additional trap 
potential that is slowly varying on the scale of the optical lattice but we neglect here 
the inhomogeneity due to such a potential and consider instead the thermodynamic 
limit. 

In terms of bosonic creation and annihilation operators, a x and a x , our Hamiltonian 
is expressed as 

H = -\ J^(4a y + a x a y ) + A ^(-l) x a x a x + U a^a^a^ - 1). (1) 

(xy) x x 

The sites x are in a cube AcZ d with opposite sides identified (i.e., a (^-dimensional 
torus) and (xy) stands for pairs of nearest neighbors. Units are chosen such that 
h 2 /m = 1. 
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The first term in ((TJ) is the discrete Laplacian ^ xy )(a x — a y )(a x — o y ) minus 
2d^ x a x a x , i.e., we have subtracted a chemical potential that equals d. 

The optical lattice gives rise to a potential A(— l) x which alternates in sign between 
the A and B sublattices of even and odd sites. The inter-atomic on-site repulsion is 
U, but we consider here only the case of a hard-core interaction, i.e., U = oo. If A = 
but U < oo we have the Bose-Hubbard model. Then all sites are equivalent and the 
lattice represents the attractive sites of the optical lattice. In our case the adjustable 
parameter is A instead of U and for large A the atoms will try to localize on the B 
sublattice. The Hamiltonian (0) conserves the particle number N and it can be shown 
that, for U = oo, the lowest energy is obtained uniquely for N = ||A|, i.e., half the 
number of lattice sites. Because of the periodic potential the unit cell in this model 
consists of two lattice sites, so that we have on average one particle per unit cell. This 
corresponds, physically, to filling factor 1 in the Bose-Hubbard model. 

For given temperature T, we consider grand-canonical thermal equilibrium states, 
described by the Gibbs density matrices Z~ x exp(—j3H) with Z the normalization 
factor (partition function) and (3 = 1/T the inverse temperature. Units are chosen so 
that Boltzmann's constant equals 1. The thermal expectation value of some observable 
O will be denoted by (O) = Z^Tr O exp( -/3H) . 

Our main results about this model can be summarized as follows: 

1. If T and A are both small, there is Bose-Einstein condensation. In this parameter 
regime the one-body density matrix 7(x, y) = (a x a y ) has exactly one large 
eigenvalue (in the thermodynamic limit), and the corresponding condensate wave 
function is 0(x) = const ant. 

2. If either T or A is big enough, then the one-body density matrix decays expo- 
nentially with the distance |x — y|, and hence there is no BEC. In particular, 
this applies to the ground state T = for A big enough, where the system is in 
a Mott insulator phase. 

3. The Mott insulator phase is characterized by a gap, i.e., a jump in the chemical 
potential. We are able to prove this, at half-filling, in the region described in 
item 2 above. More precisely, there is a cusp in the dependence of the ground 
state energy on the number of particles; adding or removing one particle costs 
a non-zero amount of energy. We also show that there is no such gap whenever 
there is BEC. 

4. The interparticle interaction is essential for items 2 and 3. Non-interacting 
bosons always display BEC for low, but positive T (depending on A, of course). 

5. For all T > and all A > the diagonal part of the one-body density matrix 
(a x a x ) (the one-particle density) is not constant. Its value on the A sublattice 
is constant, but strictly less than its constant value on the B sublattice and this 
discrepancy survives in the thermodynamic limit. In contrast, in the regime 
mentioned in item 1, the off-diagonal long-range order is constant, i.e., (a x a y ) ~ 
0(x)0(y)* for large |x — y| with 0(x) =constant. 
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Figure 1: Schematic phase diagram at half- filling 

Because of the hard-core interaction between the particles, there is at most one 
particle at each site and our Hamiltonian (with U = oo) thus acts on the Hilbert space 
7i = &) xeA C 2 . The creation and annihilation operators can be represented as 2 x 2 
matrices with 

a ^(S J)' a ^(°i o)' a ^(i o)' 

for each xG A. More precisely, these matrices act on the tensor factor associated with 
the site x while a x and a x act as the identity on the other factors in the Hilbert space 

W <8>x.A~- 

The Hamiltonian can alternatively be written in terms of the spin 1/2 operators 

The correspondence with the creation and annihilation operators is 

a x = ^x + = ) a x = ~ i^x = S x , 

and hence a x a x = + |. (This is known as the Matsubara-Matsuda correspon- 
dence J7|.) Adding a convenient constant to make the periodic potential positive, the 
Hamiltonian ([TJ for U = 00 is thus equivalent to 

H = 4E(^y +5 x 5 y + ) + A^[| + (-i) x 5g 

<xy> x 

= -E( 5 X + 5 x4) + A E[l + (- 1 ) X5 x]- (2) 

(xy) x 

Without loss of generality we may assume A > 0. This Hamiltonian is well known as 
a model for interacting spins, referred to as the XY model ^H] • The last term has the 
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interpretation of a staggered magnetic field. We note that BEC for the lattice gas is 
equivalent to off-diagonal long range order for the 1- and 2-components of the spins. 

The Hamiltonian (J2J) is clearly invariant under simultaneous rotations of all the 
spins around the 3-axis. In particle language this is the U(l) gauge symmetry asso- 
ciated with particle number conservation of the Hamiltonian (0). Off-diagonal long 
range order (or, equivalently, BEC) implies that this symmetry is spontaneously bro- 
ken in the state under consideration. It is notoriously difficult to prove such symmetry 
breaking for systems with a continuous symmetry. One of the few available techniques 
is that of reflection positivity (and the closely related property of Gaussian domina- 
tion) and fortunately it can be applied to our system. For this, however, the hard core 
and half-filling conditions are essential because they imply a particle-hole symmetry 
that is crucial for the proofs to work. Naturally, BEC is expected to occur at other 
fillings, but no one has so far found a way to prove condensation (or, equivalently, long- 
range order in an antiferromganet with continuous symmetry) without using reflection 
positivity and infrared bounds, and these require the addtional symmetry. 

Reflection positivity was first formulated by K. Osterwalder and R. Schrader 
in the context of relativistic quantum field theory. Later, J. Frohlich, B. Simon and 
T. Spencer used the concept to prove the existence of a phase transition for a classical 
spin model with a continuous symmetry (201, and E. Lieb and J. Frohlich j2Ij as well 
as F. Dyson, E. Lieb and B. Simon ^H] applied it for the analysis of quantum spin 
systems. The proof of off-diagonal long range order for the Hamiltonian (J2J) (for small 
A) given here is based on appropriate modifications of the arguments in [T5] . 

2 Reflection Positivity 

In the present context reflection positivity means the following. We divide the torus 
A into two congruent parts, A L and A R , by cutting it with a hyperplane orthogonal to 
one of the d directions. (For this we assume that the side length of A is even.) This 
induces a factorization of the Hilbert space, 7i = Hh £g> 7Yr, with 

H L , R = (g) C 2 . 

xeA LjR 

There is a natural identification between a site xG Al and its mirror image e Ar. 
If F is an operator on 7i = TC^ we define its reflection OF as an operator on T~Cr in the 
following way. If F = F x operates non-trivially only on one site, x G Al, we define 
OF = VF&xV' where V denotes the unitary particle-hole transformation or, in the spin 
language, rotation by ir around the 1-axis. This definition extends in an obvious way 
to products of operators on single sites and then, by linearity, to arbitrary operators 
on Hi,. Reflection positivity of a state ( ■ ) means that 

(FOF) > (3) 

for any F operating on Ti^. Here F is the complex conjugate of the operator F in the 
matrix representation defined above, i.e., defined by the basis where the operators 
are diagonal. 
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We now show that reflection positivity holds for any thermal equilibrium state of 
our Hamiltonian. We can write the Hamiltonian (j2J) as 



H — Hl + — \ (S^S y + Sy 

(xy)eM 



(4) 



where H L and H R act non-trivially only on 7^ L and 7Yr, respectively. Here, M denotes 
the set of bonds going from the left sublattice to the right sublattice. (Because of the 
periodic boundary condition these include the bonds that connect the right boundary 
with the left boundary.) Note that = 9Hl, and 

2J (S+S~ + S~ Sy ) = 2J (S*0S+ + S~9S~). 

(xy)GM (xy)GM 

For these properties it is essential that we included the unitary particle-hole trans- 
formation V in the definition of the reflection 9. For reflection positivity it is also 
important that all operators appearing in H (0J) have a real matrix representation. 
Moreover, the minus sign in (jlj) is essential. 
Using the Trotter product formula, we have 

Tr F9Fe~ l3H = lim Tr F9F Z„ 



with 



Z n 



&H L9e -& Bl j] ( 1 + £[s+es+ + szes;)] 



<xy)€M 

Observe that Z n is a sum of terms of the form 



(5) 



(6) 



with Ai given by either e or \J -§^S^ or \j ^-i> x . All the Ai are real matrices 



In x 



and therefore 



Tr H F9FU i Ai9A i = Tr H Fll i A i 9 FR-A,- = \Tt Hl F]JA\ > 0- 



(7) 



Hence Tr F9F Z n is a sum of non-negative terms and therefore non-negative. This 
proves our assertion. 



3 Proof of BEC for Small A and T 

The main tool in our proof of BEC are infrared bounds. More precisely, for p e A* (the 
dual lattice of A), let S* = (A)" 1 / 2 ^ x S* exp(ip • x) denote the Fourier transform of 
the spin operators. We claim that 

air. (8) 
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with E p = X)f=i(l — cos(pi)). Here, pi denotes the components of p, and ( , ) denotes 
the Duhamel two point function at temperature T, defined by 

(A, B) = f Tr (Ae- sPH Be-^ H ) ds/Ti e^ H (9) 
Jo 

for any pair of operators A and B. Because of invariance under rotations around the 
S 3 axis, (JBJ) is equally true with 5" 1 replaced by S 2 , of course. 

The crucial lemma (Gaussian domination) is the following. Define, for a complex 
valued function h on the bonds (xy) in A, 

Z(h) = Tr exp[- (3 K(h)}, (10) 

with K(h) the modified Hamiltonian 

K W = \ E ((# - S- - K y ) 2 + (S- -S;-h^) 2 )+\J2ti + (-l) x 5x] • (11) 

(xy) x 

Note that for h = 0, agrees with the Hamiltonian if, because (S^) 2 = 0. We 

claim that, for any real valued h, 

Z(h) < Z(0). (12) 

The infrared bound then follows from d 2 Z(eh)/de 2 \ £= o < 0, taking h xy = exp(ip • 
x) — exp(ip • y). This is not a real function, though, but the negativity of the (real!) 
quadratic form d 2 Z(eh)/de 2 \ e=0 for real h implies negativity also for complex- valued 
h. 

The proof of (|T2*|) is very similar to the proof of the reflection positivity property 
(£3) given above. It follows along the same lines as in ^HJ , but we repeat it here for 
convenience of the reader. 

The intuition behind (|12j) is the following. First, in maximizing Z(h) one can 
restrict to gradients, i.e., h xy = /i x — h y for some function h x on A. (This follows from 
stationarity of Z(h) at a maximizer h max ).) Reflection positivity implies that (ABB) 
defines a scalar product on operators on 7i L , and hence there is a corresponding 
Schwarz inequality. Moreover, since reflection positivity holds for reflections across 
any hyperplane, one arrives at the so-called chessboard inequality, which is simply a 
version of Schwarz's inequality for multiple reflections across different hyperplanes. 
Such a chessboard estimate implies that in order to maximize Z(h) it is best to choose 
the function h x to be constant. In the case of classical spin systems [201; ^ ms intuition 
can be turned into a complete proof of ([12)) . Because of non-commutativity of K(h) 
with K(0) = H, this is not possible in the quantum case. However, one can proceed 
by using the Trotter formula as follows. 

Let /imax be a function that maximizes Z(h) for real valued h. If there is more 
than one maximizer, we choose /i max 1° be one that vanishes on the largest number 
of bonds. We then have to show that actually /i max = 0. If /i max ^ 0, we draw a 
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hyperplane such that h Ky ^ for at least one pair (xy) crossing the plane. We can 
again write 

K{h) = K L {h) + K R {h) + \ {{S+-S y -K y y + {S^-S;-K y ) 2 ). (13) 

(xy)eM 

Using the Trotter formula, we have Z(h) = lim n ^ 00 a n , with 



OL n = Tr 



-PK L /n-f3K R /n 



n 



,-P(S+-S--h xy ) 2 /4n -{3(S--S+-h xy f/4n 



(xy)eA/ 



(14) 



For any matrix, we can write 



-D 2 



(An)- 1 / 2 [ dke [kD e- k2 /\ 
Jr 



(15) 



If we apply this to the last two factors in (|14j). and note that S y = 8S£ if (xy) G M. 
Denoting by x 1; . . . , x; the points on the left side of the bonds in M, we have that 



a n = (An)-" 1 I d Z! "kl\ 

J R 2nl 

xe -fc 2 /4 e -ifci/i xl ^ xl /3 1 / 2 /2n 1/2 



]2nl i 



PKh/n p -PK R /n ifc 1 (S+-0S+)/3 1 / 2 /2n 1 /2 



(16) 



Here we denotes k 2 = ^2 ^1 f° r short. Since matrices on the right of M commute with 
matrices on the left, and since all matrices in question are real, we see that the trace 
in the integrand above can be written as 



Tr 



-/3K L /n ifciS+^/a/anVa 



Tr 



3 -/3i^/n p ifci6>,S+/3 1 /2/2ni/2 



(17) 



Using the Schwarz inequality for the k integration, and 'undoing' the above step, we 
see that 



Kl 2 < ( {^Y nl I d 2nl ke~ k /4 

I R 2nl 



xTr \ e -PKr,M e -W K L/n e ikx(S£i -95+)/3 1 / 2 /2n 1 /2 



x ( (47i)- nl / d 2nl ke- k2/4 

;Tr \ e -P eK R/ n e -f 3K R/ n e ik i( s Zi- ss Zi)P 1/2 / 2nl/2 . 



In terms of the partition function Z(h), this means that 

\Z(h max )\ 2 < Z(hW)Z(h®), 



(19) 
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where h^ 1 ' and h^) are obtained from h ma , x by reflection across M, in the following 
way: 

{/i xy if x, y G A L 
h My ifx,yGA R (20) 
if (xy) e M 

and is given by the same expression, interchanging L and R. Therefore also 
and must be maximizers of Z(h). However, one of them will contain strictly more 
zeros than /i max , since h max does not vanish identically for bonds crossing M. This 
contradicts our assumption that h max contains the maximal number of zeros among 
all maximizers of Z(h). Hence h max = identically This completes the proof of (jl2j) . 

The next step is to transfer the upper bound on the Duhamel two point function 
(jHj) into an upper bound on the thermal expectation value. This involves convexity 
arguments and estimations of double commutators like in Section 3 in [Ej. For this 
purpose, we have to evaluate the double commutators 

[Si, [H, S^]} + [S 2 p , [H, S 2 _ p ]] =-l- ( H -lA\A\+2j2 cos p • (x - y)) . (21) 

1 1 <xy> 

Let C p denote the expectation value of this last expression, 

C P = ([^,[^5 ; -p]] + [^[^5 2 p]])>0. 

The positivity of C p can be seen from an eigenfunction-expansion of the trace. From 
[T%1 Corollary 3.2 and Theorem 3.2] and (jBJ) we infer that 




(5j5i p + S 2 p S 2 _ p ) < - coth J^C p E p /A. (22) 



Using coth x < 1 + 1/x and Schwarz's inequality, we obtain for the sum over all p ^ 0, 

YUm , + a&j < \ £ i + \ ( £ 1,2 (E c p) 1/2 - m 

p^O H p^O P p^O p p^O 

We have X]p eA , C p = —2(H) + A|A|, which can be bounded from above using the 
following lower bound on the Hamiltonian: 

H>~ Kd + l)+4A 2 ] 1/2 + |A|A|. (24) 
This inequality follows from the fact that the lowest eigenvalue of 

2d 2d 

(25) 

i=l i=l 

is given by —^[d(d + 1) + 4A 2 ] 1 / 2 . This can be shown exactly in the same way as [T%1 
Theorem C.l]. Since the Hamiltonian H can be written as a sum of terms like (|25|). 
with yi the nearest neighbors of x, we get from this fact the lower bound (J24j) . 
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With the aid of the sum rule 



peA* 



JAI 
2 



(which follows from (S 1 ) 2 = (S 2 ) 2 = 1/4), we obtain from fl2Hl and (jUj) the following 
lower bound in the thermodynamic limit: 

lim -LlS^ + S'S 2 ) 

A^oo A 

>\-\(\ [did + 1) + 4A 2 ] 1/2 c d ) V2 - ic d , (26) 



c d = 77^ I dp^r. (27) 



with Cd given by 

(27r)^y h ^ ]d -E, 

This is our final result. Note that q is finite for d > 3. Hence the right side of ()26|) is 
positive, for large enough j3, as long as 

2 1 1) 



A < -7T — 



eg 4 

In <i = 3, C3 ~ 0.505 ^H], and hence this condition is fulfilled for A < 0.960. In 
it was also shown that ded is monotone decreasing in d, which implies a similar result 
for all d > 3. 

The connection with BEC is as follows. Since H is real, also 7(x, y) is real and we 
have 

7 (x,y) = (5 x + 5 y ) = (^4 + ^4). 
Hence, if tp = |A| -1 / 2 denotes the constant function, 

(^1^) = (S^ + SlS 2 ), 

and thus the bound implies that the largest eigenvalue of 7(x, y) is bounded from 
below by the right side of (|26[). In addition one can show that the infrared bounds 
imply that there is at most one large eigenvalue (of the order |A|), and that the 
corresponding eigenvector (the 'condensate wave function') is strictly constant in the 
thermodynamic limit pQ . The constancy of the condensate wave function is surprising 
and is not expected to hold for densities different from |, where particle-hole symmetry 
is absent. In contrast to the condensate wave function the particle density shows the 
staggering of the periodic potential PjJ Thm. 3]. It also contrasts with the situation 
for zero interparticle interaction, as discussed at the end of this paper. 

In the BEC phase there is no gap for adding particles beyond half filling (in the 
thermodynamic limit): The ground state energy, E}., for ||A| + k particles satisfies 



(28) 
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(with a constant that depends on k but not on |A|.) The proof of (f2*%|) is by a variational 
calculation, with a trial state of the form (Sq ) fc |0), where |0) denotes the absolute 
ground state, i.e., the ground state for half filling. (This is the unique ground state 
of the Hamiltonian, as can be shown using reflection positivity. See Appendix A in 
PP.) Also, in the thermodynamic limit, the energy per site for a given density, e(g), 
satisfies 

e(g) - e(|) < const, (g - \f . (29) 
Thus there is no cusp at q = 1/2. To show this, one takes a trial state of the form 

IW> = e bE -^(4 + 5)|0>- (30) 

The motivation is the following: we take the ground state and first project onto a given 
direction of S 1 on some site y. If there is long-range order, this should imply that 
essentially all the spins point in this direction now. Then we rotate slightly around 
the S^-axis. The particle number should then go up by e|A|, but the energy only by 
£ 2 |A|. We refer to [H Sect. IV] for the details. 

The absence of a gap in the case of BEC is not surprising, since a gap is charac- 
teristic for a Mott insulator state. We show the occurrence of a gap, for large enough 
A, in the next section. 

4 Absence of BEC and Mott Insulator Phase 

The main results of this section are the following: If either 

• A > and T > d/(21n2), or 

• T > and A > such that A + |e(A)| > d, with e(A) = ground state energy per 
site, 

then there is exponential decay of correlations: 

7( x > y) < (const.) exp(-/c|x - y|) (31) 

with k > 0. Moreover, for T = 0, the ground state energy in a sector of fixed particle 
number N = ||A| + k, denoted by E k , satisfies 

E k + E_ k -2E >(\+\e(\)\-d)\k\. (32) 

I.e, for large enough A the chemical potential has a jump at half filling. 

The derivation of these two properties is based on a path integral representation of 
the equilibrium state at temperature T, and of the ground state which is obtained in 
the limit T — > oo. density matrix. The analysis starts from the observation that the 
density operator e~" H has non- negative matrix elements in the basis in which {S^} are 
diagonal, i.e. of states with specified particle occupation numbers. It is convenient to 
focus on the dynamics of the 'quasi-particles' which are defined so that the presence of 
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one at a site x signifies a deviation there from the occupation state which minimizes 
the potential-energy. Since the Hamiltonian is H = Hq + XW, with Hq the hopping 
term in (J2J) and W the staggered field, we define the quasi-particle number operators 
n x as: 

ala x for x even 



71* 



1 — a x a x for x odd 



Thus n x — 1 means presence of a particle if x is on the A sublattice (potential 
maximum) and absence if x is on the B sublattice (potential minimum). 

The collection of the joint eigenstates of the occupation numbers, {|{^ x })}, pro- 
vides a convenient basis for the Hilbert space. The functional integral representation 
of ({n x } \ e~ /3 ^ Ho+xw ^ \ {n x }) involves an integral over configurations of quasi-particle 
loops in a space x time for which the (imaginary) 'time' corresponds to a variable 
with period (3. The fact that the integral is over a positive measure facilitates the 
applicability of statistical-mechanics intuition and tools. One finds that the quasi- 
particles are suppressed by the potential energy, but favored by the entropy, which 
enters this picture due to the presence of the hopping term in H . At large A, the 
potential suppression causes localization: long 'quasi-particle' loops are rare, and the 
amplitude for long paths decays exponentially in the distance, both for path which 
may occur spontaneously and for paths whose presence is forced through the inser- 
tion of sources, i.e., particle creation and annihilation operators. Localization is also 
caused by high temperature, since the requirement of periodicity implies that at any 
site which participates in a loop there should be be at least two jumps during the 
short 'time' interval [0,(3) and the amplitude for even a single jump is small, of order 
(3. 

The path integral described above is obtained through the Dyson expansion 

e t(A +B) = e tA [ B(t m ) ■ ■ ■ Bitjdh ■■■dt m (34) 

m>0 Jo<h<t 2 <-<t m <t 

for any matrices A and B and t > 0, with Bit) = e~ tA Be tA . (The m = term in the 
sum is interpreted here as 1.) 

In evaluating the matrix elements of e~@ H = e _/3 ^ 0+AM/ \ in the basis {|{^ x })}, 
we note that W is diagonal and ({n x }|i?o|{w x }) is non-zero only if the configurations 
{n x } and {n' x } differ at exactly one nearest neighbor pair of sites where the change 
corresponds to either a creation of a pair of quasi-particles or the annihilation of such 
a pair. I.e., the matrix elements are zero unless n x = n x for all x except for a nearest 
neighbor pair (xy), where n x = n y , n x = n' , and n x + n x = 1. In this case, the matrix 
element equals —1/2. 
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Introducing intermediate states, the partition function can thus be written as fol- 
lows: 



with the interpretation i — t m — (3. Note that the factor in the last two lines of (|35|) 
equals (l/2) m if adjacent elements in the sequence of configurations {n^} differ by 
exactly one quasi-particle pair, otherwise it is zero. 



Expansions of this type are explained more fully in [22] • A compact way of writing 
dSU) is: 



Here the 'path' uo stands for a set of disjoint oriented loops in the 'space-time' A x [0, (3], 
with periodic boundary conditions in 'time'. Each uj is parametrized by a number of 
jumps, m, jumping times < t\ < t 2 < ■ ■ ■ < t m < (3, and a sequence of configurations 
{fix }, which is determined by the initial configuration {rix^} plus a sequence of 'rungs' 
connecting nearest neighbor sites, depicting the creation or annihilation of a pair of 
neighboring quasi-particles (see Fig. |2J). As in Feynmann's picture of QED, it is 
convenient to regard such an event as a jump of the quasi-particle, at which its time- 
orientation is also reversed. The length of u, denoted by is the sum of the vertical 
lengths of the loops. The measure v (du) is determined by (J35j) : namely, for a given 
sequence of configurations {nx'}, 1 < i < m, the integration takes places over the 
times of the jumps, with a measure (l/2) m dt 1 ■ ■ ■ dt m . 

One may note that the measure v{dui) corresponds to a Poisson process of ran- 
dom configurations of oriented 'rungs', linking neighboring sites at random times, and 
signifying either the creation or the annihilation of a pair of quasiparticles. The ma- 
trix element ({n^}\e~^ H \{n' y! .}) gets no contribution from rung configurations that are 
inconsistent, either internally or with the boundary conditions corresponding to the 
specified state vectors. A consistent configuration yields a family of non-overlapping 
loops which describe the motion of the quasi-particles in in the 'space-time' A x [0, /3). 
Each such configuration contributes with weight e~ A ' w ' to the above matrix element 
(another positive factor was absorbed in the measure v(du)). One may note that long 
paths are suppressed in the integral (J38|) at a rate which increases with A. 

Likewise, for x ^ y, we can write 




x(-i) m (K 1) }l^o|{4 m )})(K m )}|i/ |K m - 1 )}) 

x ({n^}\H \{n^}) ■ ■ ■ ({nV}\H \\{nM}) 



(35) 




(36) 




(37) 
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1 ■ ■ 1 1 ■ 1 1 1 

ABABABABA 

Figure 2: Loop gas describing paths of quasi-particles for particle number N = |A|/2 — 
1. A line on an A site means presence of a particle, while on a B site it means absence. 
The horizontal rungs correspond to hopping of a particle. 



where A^'^ denotes the set of all loops that, besides disjoint closed loops, contain one 
curve which avoids all the loops and connects x and y at time zero. The one-particle 
density matrix can thus be written 

/ (x , y) v{duj)e- x ^ 
7 x,y = A . 38 

J v(du)e A l w l 

For an upper bound, we can drop the condition in the numerator that the loops 
and the curve from x to y do not intersect. The resulting measure space is simply a 
Cartesian product of the measure space appearing in the denominator and the space 
of all curves, (, connecting x and y, both at time 0. Denoting the latter by S(x, y), 
we thus get the upper bound 

7 (x,y)< / v(dC)e- x ^. (39) 

The integral over paths is convergent if either A or T is small enough, and away from 
the convergence threshold the resulting amplitude decays exponentially. A natural 
random walk estimate, see PQ Lemma 4], leads to the claimed exponential bound 
provided 

d (1 - e~ px ) < A. (40) 

This includes, in particular, the cases T > d for any A, and A > d for any T. 
Exponential decay actually holds for the larger range of parameters where 

d (1 - e-^-V) < A - /, (41) 

where / = f(/3, A) = — ((3\ A|) _1 In Tr e~^ H is the free energy per site. Note that / < 0. 
This condition can be obtained by a more elaborate estimate than the one used in 
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obtaining (jHOj) from (|38j) . as shown in PJ Lemma 3]. The argument there uses reflection 
positivity of the measure v(du). Using simple bounds on / one can then obtain from 
(j4*Tj) the conditions stated in the beginning of this section. 

The proof of the energy gap is based on an estimate for the ratio ^^I-ph where 
Vk projects onto states in Fock space with particle number N — ^\A\ + k, expressing 
numerator and denominator in terms of path integrals. The integral for the numerator 
is over configurations uj with a non-trivial winding number k. Each such configuration 
includes a collection of 'non-contractible' loops with total length at least f3\k\. An 
estimate of the relative weight of such loops yields the bound 



< (const. )(|A|/|A;|) |A;| (e 1 ^ 0118 (42) 



TrV e-' 3H 
which gives for f3 — > oo 

E k - E > (const. )|Jfe| (43) 
independently of |A|. We refer to [T] for details. 



5 The Non-Interacting Gas 

The interparticle interaction is essential for the existence of a Mott insulator phase for 
large A. In case of absence of the hard-core interaction, there is BEC for any density 
and any A at low enough temperature (for d > 3). To see this, we have to calculate the 
spectrum of the one-particle Hamiltonian — |A + V(x), where A denotes the discrete 
Laplacian and V^(x) = A(— l) x . The spectrum can be easily obtained by noting that 
V anticommutes with the off-diagonal part of the Laplacian, i.e., {V, A + 2d} = 0. 
Hence 

(-±A -d+ K(x)) 2 = (-±A - d) 2 + A 2 , (44) 
so the spectrum is given by 

d± vfecos^f + A 2 , (45) 

where p G A*. In particular, E(p) — E(0) ~ \d(d 2 + A 2 ) _1 / 2 |p| 2 for small |p|, and 
hence there is BEC for low enough temperature. Note that the condensate wave 
function is of course not constant in this case, but rather given by the eigenfunction 
corresponding to the lowest eigenvalue of — |A + A(— l) x . 



6 Conclusion 

In this paper a lattice model is studied, which is similar to the usual Bose-Hubbard 
model and which describes the transition between Bose-Einstein condensation and 
a Mott insulator state as the strength A of an optical lattice potential is increased. 
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While the model is not soluble in the usual sense, it is possible to prove rigorously all 
the essential features that are observed experimentally. These include the existence of 
BEC for small A and its suppression for large A, which is a localization phenomenon 
depending heavily on the fact that the Bose particles interact with each other. The 
Mott insulator regime is characterized by a gap in the chemical potential, which does 
not exist in the BEC phase and for which the interaction is also essential. It is possible 
to derive bounds on the critical A as a function of temperature. 
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